Constrained deformation of a confined solid: a strain induced crystal-smectic 

transition 
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We report results of computer simulations of two-dimensional hard disks confined within a quasi 
one-dimensional "hard-wall" channel, a few atomic radii wide. Starting from a commensurate tri- 
angular solid a rescaling of the system size parallel to the channel length introduces a rectangular 
distortion of the solid which, beyond a critical limit, phase separates into alternating bands of solid 
and smectic phases. The resulting solid- smectic interfaces are broad and incorporate misfit dislo- 
cations. The stress-strain curve shows large plastic deformation accompanying the crystal-smectic 
transition which is reversible. The smectic phase eventually melts into a modulated liquid with a 
divergent Lindemann parameter. 



Studies of small assemblages of molecules with one or 
more dimensions comparable to a few atomic spacings are 
significant in the context of nano-technology 1]. Design- 
ing nano-sized machines requires a knowledge of the me- 
chanical behavior of systems up to atomic scales, where, 
a priori, there is no reason for continuum elasticity theory 
to be valid. In most cases, however, the effects of finite 
size are relatively mild, showing up mainly as a vari- 
ation of the numerical value of the elastic constants 
as a function of length scale. In this Letter, we show, 
on the other hand, that small size and hard constraints 
can produce essentially new phenomena without a coun- 
terpart in the bulk system. We perform computer sim- 
ulations of the simplest possible, nontrivial, molecular 
system, namely, two-dimensional hard disk "atoms" con- 
fined within a quasi one-dimensional channel; the physics 
of which is entirely governed by geometry. Bulk hard 
disks in two dimensions are known to melt 0, 0, P| from 
a high density triangular lattice to an isotropic liquid 
with a narrow intervening hexatic phased, 0]. In con- 
trast, for channel widths of a few atomic spacings, we find 
evidence for a smectic phase which nucleates as promi- 
nent bands within the solid. The smectic phase arises 
when the size of the system in the direction parallel to 
the fixed walls is increased. A crystal to smectic tran- 
sition, though predicted for anisotropic molecules is 
unusual for hard disks - the anisotropy in this case aris- 
ing purely from the external confining potential. In this 
respect, our results resemble the phenomenon of laser in- 
duced freezing^ where an external modulated electric 
field produced by crossed laser beams induces a series of 
phase transitions 0,01 involving triangular solid, modu- 
lated liquid as well as smectic phases. The analog of the 
modulating potential in our case is, of course, provided 
by the walls which induce a periodic potential of mean 
forceiniEl 

decaying with distance from the walls. The 
nature of the crystal- smectic transition in our system, as 
we show below, are, however, different. Our results may 
be directly verified in experiments on sterically stabilized 
"hard sphere" colloids [l3| confined in glass channels and 
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FIG. 1: Results of NVT ensemble Monte Carlo simulations 
of A*' = rix X riy = 65 X 10 hard disks confined between 
two parallel hard walls separated by a distance Ly = 9.001 d 
where d is the hard disk diameter. Inset shows the geome- 
try used; only four layers have been shown for clarity. The 
reciprocal lattice vectors (RLVs) Gi and G2, the rectangu- 
lar unit cell and the lattice parameters and ay are also 
shown. At ?7 = 0.85 we have a strain free triangular lat- 
tice. Plots show pGi,i = l(-|-),2(o) the structure factor for 
RLVs Gi{ri), averaged over symmetry related directions, as 
a function of r;. The value of ri is changed by changing the 
length of the box Lx while keeping Ly and N fixed. For each 
T], the system was equilibrated over 10® Monte Carlo steps 
(MCS) and data averaged over a further 10® MCS. pci is 
non-zero throughout implying long ranged orientational or- 
der. In contrast, pa^ jumps to zero at 77 = t^ci ~ .77. Also 
plotted in the same graph is the Lindemann parameter l{0) 
which diverges below 77 — r^c^ — .7 < ri^. Note that at 77^ 
X ~ 9.54 ~ ni{— 10) — 1/2. The lines in the figure are a guide 
to the eye. 



may also be relevant for similarly confined atomic sys- 
tems interacting with more complex potentials. 

The bulk system of hard disks where particles i and j, 
in two dimensions, interact with the potential Vij — 

< d, where d is the 
r; the relative posi- 
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FIG. 2: A plot of the normal stress versus the conjugate 
strain Ed ~ {rjo — v)/'no ivo ~ 0.85) obtained from our Monte 
Carlo simulations of 65 x 10 hard disks showing a typical Van 
der Waals loop in the constant strain ensemble. Data for 
the plot is obtained by equilibrating at each strain value for 
2 X 10* MCS and averaging the data for a further 3 x lO" MCS. 
The stress for the hard disk system has been calculated by the 
standard method fl5|l by averaging the collision probability. It 
is minimum at 77 = 77C2 ~ -74. The entire cycle consisting of 
increasing £d(o) and again decreasing to zero (+) using typi- 
cal parameters appropriate for an atomic system, corresponds 
to a real frequency of ~ lOOKHz. The lines in the figure are a 
guide to eye. The solid (ry > rjc^ ) , two phase (r^ca < Tj < rjc^ ) 
and smectic (Sm) (t; < r^ca) regions are indicated in the fig- 
ure. We have repeated this calculation with a cycle frequency 
lOK Hz — IM Hz with no essential change in the results. 



tion vector of the particles, has been extensively^ 0, Q 
studied. Apart from being easily accessible to theo- 
retical treatment I experimental systems with nearly 
"hard" interactions viz. sterically stabilized colloids|l3j 
are available. The hard disk free energy is entirely en- 
tropic in origin and the only thermodynamically relevant 
variable is the number density p — N/V or the packing 
fraction 77 = (7r/4)/5d^. Simulations 4] , experimental'lSj 
and theoretical 14] studies of hard disks show that for 7/ > 
77^ = .719 the system exists as a triangular lattice which 
transforms to a liquid below rjm = .706. The small inter- 
vening region contains a hexatic phase predicted by the 
Kosterlitz-Thouless-Halperin-Nelson- Young theory 5] of 
two dimensional melting. The surface free energy of the 
hard disk system in contact with a hard wall has also 
been obtained [T3| taking care that the dimensions of the 
system are compatible with a perfect (strain-free) trian- 
gular lattice. The effect of in-commensuration, though, 
which may routinely arise in experiments, surprisingly, 
have not received the attention it deserves. We attempt 
to address this aspect as follows. 

Consider a narrow channel in two dimensions of width 
Ly defined by hard walls at 77 = and Ly (Vivaii(r) = 



for < r„ 



< Ly and 



00 otherwise) and length with 



the direction x implying x + = x. In order that the 
channel may accommodate ni layers of a homogeneous, 
triangular lattice with lattice parameter oq of hard disks 
of diameter d, (Fig^l one needs. 



Ly = — (ti; - l)ao + d 



(1) 



Defining x = 1 + 2{Ly — d)/\/3ao, the above condition 
reads x = integer = 77; and violation of Eqn.ljTjl implies 
a rectangular strain away from the reference triangular 
lattice of ni layers. The lattice parameters of a centered 
rectangular (CR) unit cell are Qx and Uy (Fig. ^ inset). 
In general, for a CR lattice with given Ly we have, ay — 
2{Ly — d)/{ni — 1) and, ignoring vacancies, — 2/ pay. 
The normal strain Sd = £xx — £11 



£d = 



ni-l 



zyy is then, 
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where the number of layers n; is the nearest integer to x 
so that Ed has a discontinuity at half -integral values of 
X- For large Ly this discontinuity and Ed itself vanishes 
as 1/Ly for all 77. 

We study the effects of this strain Ed on the hard disk tri- 
angular solid at fixed Ly large enough to accommodate a 
small number of layers ni ~ 9 — 25. The strain Ed is im- 
posed by expanding the dimension of the system L^ par- 
allel to the walls keeping Ly fixed so that Ed — (va^v) I Voj 
where rjo is the packing fraction corresponding to an un- 
strained triangular solid. We monitor the Lindemann pa- 
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Lx ^ Ly. Periodic boundary conditions are assumed in 



where the angular brackets denote averages over config- 
urations, i and j are nearest neighbors and u^i is the 
a-th component of the displacement of particle i from 
it's mean position. The parameter I diverges at the melt- 
ing transition . We also measure the structure factor 

PG = (172 E5=i®^p(~**^-''»j)) ' ^^"^ ^ = ±G!-i(77), 
the reciprocal lattice vector (RLV) corresponding to the 
set of close-packed lattice planes of the CR lattice per- 
pendicular to the wall, and ±62(77) the four equivalent 
RLVs for close-packed planes at an angle { — n/3 and 
2tt/3 in the triangular lattice) to the wall (see Fig. 
inset). 

We find, throughout, pG-> < PGi 7^ 0, a consequence of 
the hard wall constraint 12] which manifests as an oblate 
anisotropy of the local density peaks in the solid. As 77 is 
decreased (see Fig. n]for details) both and pc^ show 
a jump at 77 = 77^ close to x ~ ^ 1/2- For 77 < 77^ 
we get = with pci 7^ signifying a transition 
from crystalline to smectic like order. The Lindemann 
parameter I remains zero and shows a divergence only 
below 77 = T]c3i~ Vm) indicating a finite-size- broadened 
melting of the smectic to a modulated liquid phase. We 
have also calculated the normal stress ad = cTxx — o'yy 
(see Fig. For 77 = 770 the stress is purely hydrostatic 
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with cr^cj; = ayy as expected. As 77 decreases, the stress 
increases Unearly in the elastic limit, flattening out at 
the onset of non-linear behavior at 77 rjc^ . At 77^ , Od 
decreases and eventually becomes negative. On further 
decrease in 77 below 77C2 (Fig. <Jd approaches from 
below thus forming a Van der Waals loop typical of the 
constant strain ensemble. If the strain is reversed by 
increasing 77 back to 770 the entire stress-strain curve is 
traced back with no remnant stress at 77 = 770 showing 
that the plastic region is reversible. As Ly is increased, 

77ci merges with 77C3 for x ~ 25. If instead, and Ly 
are both rescaled to keep x fixed or periodic boundary 
conditions are imposed in both x and y directions, the 
transitions in the various quantities occur approximately 
simultaneously as expected in the bulk system. Varying 
rix in the range 10 — 1000 produces no essential change 
in results. 

For rjc3 < 1] < rjc^ we observe that the smectic order 
appears within narrow bands (Fig. |3J| which nucleate 
at 7/ci . Inside these bands the number of layers is less 
by one and the system in this range of 77 is in a mixed 
phase. A plot (Figl^^a)) of xi^)^ obtained by averaging 
the instantaneous ay from particle configurations over a 
strip spanning Ly, shows bands in which x is less by 
one compared to the crystalline regions. Once nucleated 
narrow bands coalesce to form wider bands, the dynamics 
of which is, however, extremely slow. The bands grow 
as 7] is decreased. Calculated diffraction patterns (Fig. 13 
(a)) show that, locally, within a smectic band 2> pc^ 
in contrast to the solid region where ~ PG2 7^ 0- 

Strong finite size corrections makes a complete theoreti- 
cal treatment of this problem difficult. However, signifi- 
cant progress may be made using qualitative arguments 
as we show below. The total free energy of the system 
may be decomposed as, 

T^{rj,x)^K''{v)eUx)+^''{v) (3) 

where K^{rf) is an elastic constant and T'^{rf) the free 
energy of the perfect triangular lattice in contact with a 
hard wall^3| ^tt packing fraction 77. It is clear that JF^ 
has minima for all x = • For half integral values of x 
the crystalline structure is metastable with respect to an 
intervening smectic when adjacent local density peaks of 
the solid overlap in the x direction. This overlap is fa- 
cilitated by the (oblate) anisotropy of the solid density 
peaks. For large Ly the minima in T"^ merge to pro- 
duce a smooth free energy surface independent of x- For 
small Ly all regions of the parameter space correspond- 
ing to non-integral x a-^e also globally unstable as belied 
by the loop in the stress-strain curve (Fig[21). The sys- 
tem should therefore break up into regions with various 
ni. Such fluctuations are, however, suppressed due to the 
structure of interfaces between regions of differing ni . A 
superposition of many particle positions near such an in- 
terface (see Fig. Olb)) shows that: (1) the width of the 
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FIG. 3: (color on-line) (a) Plot of x{^) ^.s a function of the 
a;/d at 77 = .76 after 3 x 10^ 5 x lO'^ and 2 x lO'' (top) MCS 
for N = 10*. Note that x = 10 in the solid and = 9 in the 
smectic regions. Arrows show the coalescence of two bands as 
a function of time. The panel on the right shows calculated 
diffraction patterns for the solid (top), smectic (middle) and 
inter-facial (bottom) regions. Inversion symmetry is broken 
at the interface because of the inter-facial dislocation (see 
below) . (b) Close up view of a crystal-smectic interface from 
superimposed positions of 10^ configurations at 77 = .77. The 
colors code the local density of points from red (high) to blue 
(low). Note the misfit dislocation in the inter-facial region. 

interface is large, spanning about 10 — 15 atomic spac- 
ings and (2) the interface between ni layered crystal and 
71,; — 1 layered smectic contains a dislocation with Burger's 
vector in the y direction which makes up for the differ- 
ence in the number of layers. Note that the presence of 
these dislocations breaks inversion symmetry as observed 
in the local diffraction patterns. The core region of this 
dislocation extends over many layers with some of the 
disks within the interface alternating between positions 
corresponding to either a smectic or a solid. Each band 
of width s is therefore held in place by a dislocation- 
anti-dislocation pair (Fig. O. In analogy with classical 
nucleation theory^^, the free energy Fi, of a single band 
can be written as 

Ffa = -AFs + Ec + ^h^K^ log — (4) 
47r ao 

where h = ay/2 is the Burger's vector, AF the free en- 
ergy difference between the crystal and the smectic per 
unit length and Ec the core energy for a dislocation 
pair. Bands nucleate when dislocation pairs separated 
by s > j^b'^K^/AF arise due to random fluctuations. 
Band coalescence occurs by diffusion aided dislocation 
"climb" which is extremely improbable in a high den- 
sity phase leading to slow kinetics. The growing smectic 



4 



bands are in a state of tension in the y direction imply- 
ing Gyy > Gxx which is countered by the compressive 
stress in the crystalhne region. Growth of smectic bands 
therefore reduces Ud which attains a minimum at rj = rjc^ 
when a single band spans the entire length. Subsequently 
ad ~> 0, the value in the liquid phase. Since orienta- 
tion relationships between the crystal and smectic are 
preserved, the stress-strain relationship is completely de- 
termined by the amount of the co-existing phases which 
explains the reversibility^^. For large values of Ly the 
smectic phase vanishes since the strains involved them- 
selves go to zero. Nevertheless, transition between ni 
and ± 1 layered crystals have been observed by us. We 
must mention here that the choice of the ensemble, viz. 
constant strain, is crucial since, in the constant stress en- 
semble, the hard disk system fails at 77^ and the strain 
diverges producing a homogeneous low density gas with 
no interface. Finally, our smectic bands are reminiscent 
of "slip" or "deformation" bands which arise during plas- 
tic flow of macroscopic ductile materials ,2Ql|. 

Apart from constrained hard sphere colloids where 
our results are directly testable, strain induced crystal- 
smectic transition may be observable in experiments on 
the deformation of mono-layer atomic nano-beams or 
strips of real materials confined to lie within a channelP|. 
This is because the constraints we choose to study are 
geometrical and would exist in any system. The oc- 
currence of a smectic phase may be important for the 
tribologicalf^H properties of nano-scale machine parts 
where atomic friction plays an important role in their 
function. In the future we would like to study the ki- 
netics of the crystal- smectic transition in detail as well 
as the effect of substrate disorder on the nature of this 
transition. 
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